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Abstract The paper discusses some interesting, mainly 
philosophical paradigms of the modeling and control areas, 
which are still partly unsolved and/or only partially studied. 
Lirst the possible introduction of a prejudice free control - 
similar to the term for the modeling introduced by Rudi 
Kalman - is investigated. Next the real constraints in real 
control systems are discussed. It seems that these are the 
amplitude limit for the actuators in practical systems. Then 
the application of the HElSENBERG-type uncertainty 
relationship in control is discussed, combining the 
robustness of control and quality of modeling. Linally a 
special irregularity in the classical LQR control 
methodology is treated investing some unreachable poles 
in the anomaly. The paper constitutes a review of a few 
specific problems in the control theory. They are 
enumerated in keywords and discussed on the basis of the 
references. The authors aim at least to invite further authors 
to continue such a discussion. 

Keywords Prejudice Free Control, Irregularities in 
Classical Control Methods, Real Constraints in Control 
Systems, the Heisenberg Type Uncertainty Relation of 
Control 


1. Introduction 

Some scientists believe that everything has been solved 
in control, consequently nothing remained to study and/or 
research. The purpose of this paper is to recall some 
interesting philosophical paradigms in the areas of 
modeling and control to prove the contrary. 

Only a few questions are discussed here, but there are 
many. Our aim is to encourage scientists to find further 
unsolved problems, blazes and interesting paradigms partly 
based on the modeling and control literature, partly on 
other disciplines. 

If we can invite only a few further authors to continue 
our discussions then this effort is worth while. 

In the sequel the YOULA parameterization [1], [2], [4], [5] 


will be used to discuss regulator and control system design. 
We found that this is very good basis for education, too. 
The YOULA Parameterization 

The YOULA- (Tor Q) -parameterization is a classical 
method for linear time invariant control system to 
characterize all realizable stabilizing regulators ( ARS) by 


C = 


Q 

l-QP 


(l) 


for open-loop stable plant P e S , where S is the closed 
set of all stable proper real-rational systems, having all 
poles within the closed unit disc. The parameter” 


Q = 


c 

1 + C P 


Q eS 


( 2 ) 


ranges over all proper (Q (co = oo) is finite), stable transfer 
functions [1], [5]. Observe that Q is the transfer function 
from the r to u in the closed-loop (see Fig. 1), where 
y n is the output disturbance (or noise) signal in a SI SO 
(Single /nput Single Output) system. 



Figure 1 . Closed-loop with an ARS regulator 

The transfer characteristics of the closed-loop can be 
easily computed 


y = QPr-(l-QP)y n = y t + y d (3) 

where y t is the tracking (servo) and y d is the regulating 
(or disturbance rejection) independent behaviors of the 
closed-loop response, respectively. 

Because the ARS regulator represented in Fig. 3 was 
formulated for an one-degree of freedom ( IDF) control 
system, it is not surprising that the tracking part y t of the 
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transfer characteristics between y and r can not be set 
independently of the regulating behavior y d , i.e. 
independently of Q . 



Figure 2. The modified control system with an ARS regulator opening 
the closed-loop 

The Y-parameterization ’’almost” opens the closed-loop. 
Here ’’almost” means that yt=QPr is obtained instead 

of a real open-loop case with y t = Pr . So we need a Q~ l 
prefilter shown in Fig. 2, when the ARS regulator really 
’’virtually” opens the closed-loop as 

y = Pr-(\-QP)y n =y t +y d (4) 



Figure 3. The K-B-parameterized 2DF system with an ARS regulator 

An important and new observation of the authors was 
that the scheme in Fig. 2 is equivalent to the special control 
system given in Fig. 4 and its parameterization has been 
named as Keviczky-B any asz-(KB) parameterization [1], [2]. 
Since in the case of the special structure presented in Fig. 3 
we have y t =Pr , i.e., (3) holds, it is easy to introduce a 
new general form of any 2DF control systems providing 

y = QvPyr~(}-QP)yn=yt+yd (5) 

if a serial compensator Q r is applied additionally as the 
Fig. 5 shows. 



Figure 4. The general form of the K-B-parameterized 2DF control 
system 

Here and in the sequel the general notation y r will be 

used for the reference signal for general 2DF systems. 
Equation (5) shows that the tracking properties 
yt = Qr Pyx can independently be designed from the 
regulating behavior y d = (i -QP)y a by Qr- 

The last scheme was later named as a generic 
two-degree of freedom ( G2DF) system [1], [2]. The K-B 
parameterization for closed-loop control is not so widely 
known as the Youla-Kucera- ( Y-K) parameterization [4] 


however, it is much closer to a control engineering view 
and its most important advantage in 2DF systems is that it 
virtually opens the closed-loop. However, this 
parameterization can only be applied for open-loop stable 
processes. 



Figure 5. The generic 2DF ( G2DF) control system 

A G2DF control system is shown in Fig. 5, where 
y r ,u,y and y n are the reference, process input, output 
and disturbance signals, respectively. The optimal 

discrete-time ARS regulator of the G2DF scheme [1], [2] is 
given by 

C R " K " g ° R " G " P; ' (6) 

° \-RK n P 1 -QP l -R n G n P_z~ d 

where 

Qo=Q n = R n r n = R n G n p ; 1 (7) 

is the associated Y-parameter [1], furthermore 
Q t = R r K r = R r G r P+ l 

K n =G n P;' ; K T =G t F? (8) 

assuming that the continuous-time process is factorable as 

p=p_p_ = p_py sT ' (9) 

and a discrete-time process is factorable as 

G = G + G_ = G + G_z~ d (10) 

where P + , G + means the inverse stable ( IS) and P -, G_ 
the inverse unstable (JU) factors, respectively. Here T d is 

the continuous time delay and z~ d corresponds to the 
discrete time delay, which is the integer multiple of the 
sampling time T s . 

It was shown [1], [2] that the optimization of the G2DF 
scheme can be performed in H 2 and norm spaces by 
the proper selection of the serial K r and embedded K n 
filters (compensators). These optimizations are reduced to 
the optimal computation of the G v and G n embedded 
filters. If G r and G n are optimally selected, then C Q 
denotes the optimal ARS regulator in (6). It is interesting to 
see how the transfer characteristics of this system look like: 

y = R r K r Py r -(\-R n K n P)y n = 

r / t \ ( 11 ) 

= R r G I P_e~ sT *y I -(l-R a G n P_e~ s fy n = y t +y d 
or 
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y = R I K I Gy r -(l-R n K B G)y B = 

, ( a\ ( 12 ) 

= R T G r G_z- d y I -(l-R B G n G_z- d )y n =y t +y d 

Here R r and are stable and proper transfer 
functions, that are partly capable to place desired poles in 
the servo and the regulatory transfer functions, furthermore 
they are usually referred as reference signal and output 
disturbance predictors. They can even be called as 
reference models, so reasonably R r (oo = 0) = 1 and 
R n (oo = 0) = 1 are selected. In this case the obtained 
regulator is always an integrating one. 


2. Prejudice Free Control 


The knowledge of a process is never exact, independent 
of the method how its model is determined, whether 
measurement-based identification (ID) or 
physico-chemical theoretical considerations are used. The 
uncertainty of the plant can be expressed by the absolute 
model error 


AP= P-P 


and the relative model error 


l 


AP 

P 


P-P 

P 


(13) 


(14) 


where P is the available nominal model used for 
regulator design and P is the real plant. 

The parameters of the plant may change in terms of their 
nominal values in a given range. The closed-loop control 
system needs to be stable under the given uncertainty 
ranges of the parameters. 

Suppose that the open loop is stable. The regulator 
designed for the nominal plant ensures the stability of the 
nominal closed-loop control system. Let us analyze 
whether the system remains stable with the parameter 
uncertainties of the open loop. Stability is maintained if the 
NYQUIST diagram of the modified open loop does not 
encircle the -1 + 0/ point. 

If there is an uncertainty AP (or parameter change) in 
the transfer function of the plant, then if we apply the same 
regulator this uncertainty appears in the absolute error 
A L = CAP of the loop transfer function, whereas its 
relative model error is 


AL _ L-L _ CP-CP _ P-P 
L L CP P 


(15) 


Here L denotes the nominal and L denotes the real 
loop transfer function. 



Figure 6. Change in the NYQUIST diagram of an uncertain system 


Robust stability means that the closed-loop control 
system should not display unstable behavior even in the 
“worst case” parameter changes. The bound for AL can 
be formulated based in Fig. 6 by taking the simple 
geometrical considerations into account: the NYQUIST 
diagram will not encircle the -l + 0y point, if the 
following relationship is satisfied for all frequencies: 

|a£,(y'(o)| = |^(yco)| |Z(yco)| < |i + Z(_/co)| Vco (16) 

With farther straightforward manipulations the 
necessary and sufficient condition for robust stability is 
obtained as 


or 


£(j(o) < 


1 + L(ja>) 

UH 


Hj*) 


Voo 


f(j(o) £(j(o) <1 Voo 


(17) 


(18) 


where T = Lj{\ + L j is the nominal complementary 
sensitivity function. 

This form is also called the dialectic relationship of 
robust stability. In the design the first factor |T (yoo)| is 

calculated for the supposed (known) nominal parameters of 
the plant, and thus it depends on the designer. The second 

factor does not (or only partly) depends on the designer, 

as it contains the uncertainties in the knowledge of the 
plant or its unexpected parameter changes. In those 
frequency ranges where the uncertainty is large, 
unfortunately only small transfer gain can be designed for 

the closed-loop. Where |r(yco)| is high, very accurate 

information is needed to reach a small error. The higher the 
absolute value of the complementary sensitivity function, 
the smaller the permissible parameter uncertainty. 
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The robust stability condition (4), (5) and (6) can be 
rearrange in the form of 



i.e., it does not depend on the model P but only on the 
reference model or the design goal. 

The reference model is an important parameter of the 
general YOULA design, by means of which the condition of 
robust stability (22) can be guaranteed. Let 


where the function is plotted in Fig. 7. The 

interpretation of this function is very interesting. For small 
modeling error a model based controller design is 
suggested, which usually based on the inverse of the 

nominal model P~ l . For very large errors no regulator 
design is advised. However, in the midrange domain, 
where the error is around 100 %, the regulator design 
practically does not depend on our knowledge of the 
process. This area can be called “prejudice free ” domain 



R 


n 


1 +sT 


(23) 


then the constraining condition of the right side of (22) can 
be seen in Fig. 8. Given the latter condition and choosing 
first-order reference model R n , we see that robust stability 

can be ensured even in the case of 100 % relative model 
error. Furthermore for the high frequency domain a real 
prejudice free case is obtained, 



Figure 7. The r\ in the function of l 


Figure 8. Condition constraining the relative model error in the case of 
the first-order reference model 


Prejudice Free Control for YOULA Parameterized 
Systems 

The condition of robust stability for the YP control loops 
can be further simplified so the expression (18) becomes 


\QPt\ = \R n G n p;'pt\ = 

= |i? n G n A^| = kc? n / , _|w<i 


R n Gj_e' T H 


Vco 


( 20 ) 


where T d is the dead time of the model and 



The inequality (17), limiting the relative error, is now 


d-/®) < 



Vco 


( 21 ) 


If the process is IS, i.e., P_ = 1, then G = 1 can be 

chosen and the condition of robust stability can be further 
simplified as 


I'O'AH^T Vco (22) 

n 1 n| 


If the process is IU, even the factor G n iL appears in 

(21), can significantly modify (22). Fig. 9 shows the case 
when two unstable zeros seriously decrease the prejudice 
free character of the stability. The worst case is when this 
factor has a large value in the region of the cut-off 
frequency. 



Figure 9. Conditions constraining the relative model error in the case of 
two unstable zeros 

Kalman was who tried to investigate the possibility of a 
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prejudice free identification/modeling methodology [6]. 
He could not find any general applicable results, however 
many interesting, almost philosophical statements were 
developed. 


3. What are the Real Constraints in 
Control Systems 

Process Behaviors Constraining the Control 

It was always an important question in the control theory 
papers what are the real constraints strongly effecting the 
result of the controller design. 

We will discuss here the invariant process factors and 
the actuator operating signal domain here below. These 
factors are independent of the designer and of the available 
methods being either theoretically or experimentally 
founded. 

One class of the constraints, which appear in the 
controller design is that the invariant factors of the real 
process, i.e., the unstable zeros and the time delay can not 
be eliminated with any control algorithm. So the best 
reachable closed-loop performance partly depends only on 
the process itself. If someone wants to change these 
elements only the redesign and rebuilt of the technology 
helps. 

The ideal design goals formulated by the reference 
models R and R can be reached only for stable and 
inverse stable delay free processes. In case of inverse 

-sT 

unstable processes only the approximate R r G v P e d and 

—sT 

R G P_e d goals can be reached [1]. For discrete time 

systems these goals are R r G r G z~ d and R n G n G_z~ d , 
respectively. 

The influence of the unstable zeros can be somewhat 
attenuated using the G r and G n embedded filters. So 

only in case of stable process zeros we can obtain optimal 
controller independent of invariant process factor(s). 

So before designing a control systems it is better to 
clarify what the process will allow us to reach and the final 
result will not depend on our skill or on the applied 
theoretical approach and/or methodology. 

Actuator Behaviors Constraining the Control 

There is another class of constraints which is 
independent on the applied regulator design and/or control 
systems applied. These are the always existing amplitude 
limits in the real actuators. Sometimes the theoretical 
control people forget about these constraints, because in 
their platforms, which is almost always computer 
technology, their life is limited only by the very big 
floating point number representations available at different 
software environments. Unfortunately the real actuators do 
not know these internal representation in spite of the fact 


that the modern control algorithms are based on discrete 
time computer control nowadays. The amplitude limit for a 
real actuator is always finite and we must not forget this 
reality. If we want to speed up a slow process by modern 
control, the result does not depends on the theoretical 
strength of the applied algorithm, it depend primarily on 
the available amplitude limit which can be applied at the 
output of the actuator. 

The regulators always invoke zeros to accelerate the 
process. To demonstrate this let us investigate the case 
shown in Fig. 10, where a phase-lead element is connected 
in serial before the first-order lag element. In the first 
moment, a signal value of 10 appears at the output of the 
phase-lead element and at the input of the first-order lag 
element for the effect of a unit step signal. The first-order 
lag starts with a high gradient in order to reach this value as 
soon as possible with its time constant and by the time the 
input signal is settling down the output has almost reached 
its steady-state value. The cost of the acceleration is the 
so-called over-excitation, i.e., the ratio of the initial and 
final signal values at the input of the lag. The acceleration 
can be reached only by over-excitation greater than one. In 
many cases, it is worth applying pole cancellation, when 
zeros are invoked to cancel the undesirable poles which 
cause slow operation, and instead a pole ensuring more 
favorable behavior is inserted. 



Figure 10. Insertion of a zero may accelerate the system at the cost of 
over-excitation 


Thus it is obvious that the over-excitation means the 
control equipments will have an initial peak at their output 
as a response to unit step commands or disturbances. The 
problem is that the output of the regulators in the 
closed-loop control, or the output of the actuators gaining 
the signal for the proper level are always 
amplitude-restricted. 

K/fVnax (24) 


The sensitivity function of the real closed-loop can be 
written in the following decomposed form: 


J perf 


5 = (yA) + K - t)- (t-- t) = 5 des + S real + S mod = 

^ es *~Veal *^mod 

=+(« n -n=4 s +v=i | - r ) + $ 

^des 

= M 


-'mod 


I ^mod ^contr ^mod 


(25) 


S(\eK +^re 
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Here S des = (l - R n ) is the design loss, 
^real = K -jj is the realizability loss, 
S mod - ~(t-T^-T-T refers to the modeling loss in the 

sensitivity function. In its other form S contr =^l-rj 
means the decomposed term referring to the control loss, 
^perf = (^n T )to the performance loss. Each term can be 

interpreted and explained very easily. 

The above triple decomposition of the sensitivity 
functions gives a good insight into the limit-optimality 
(limits of the optimality) of closed-loop control systems, 
i.e., the characterization of the best control achievable. As 
regards this distinction optimality criteria need to be 
created for each term, i.e., 


4st^dV^al + -4,d=ll* 


J <J X +J V +J X 

J track J des ' J real ' J mod 


\K\ 

+ 

Keall 

+ | 

l^mod | 

| Ties 

+ 

^real 

+ 

C r 

°mod 


(26) 


both for the tracking and disturbance rejection behaviors. 
Here the notation |^]|| is used to express the optimality 

criterion. In strict mathematical analysis this notation is 
used to refer to the chosen norm of the transfer function. It 
is not an easy task to optimize all three terms 
simultaneously. In practice iterative techniques are used, 
whereby the optimization problem is solved step-by-step. 

The optimization of the first term in the decomposition 
of the sensitivity function (25) means the determination of 

the best (fastest) reachable reference models R x = A° pt 

and R n = R° pt , i.e., the solution of the task under the 
constraints 


simple cases. The optimal reference models are usually 
determined by simulation (CAD tools). Note that under the 
given constraints faster reference models cannot be used to 
solve the task (27). Vice versa, if under the constraints and 
design goal there is no solution for the reference models 
then the only option is to prescribe a less demanding 
(usually slower) model. Thus the best (fastest) reachable 
response of the closed system basically depends on the 
constraints of the control output. Of course, equation (27) 
contains the applied regulator and also the process in a 
complex way; thus it is a closed-loop. Therefore the 
optimality of the regulator depends on the process, the 
model and the invariant factors. 

Even in the case of a very careful design it can happen 
that the over-actuated output of the obtained regulator is 
beyond the acceptable signal domain. Then the original 
design goals need to be reduced. The advantage of the KB 
parameterization of generic 2DF control loops is that only 
the problematic (over-demanding) reference models R r 

and R n need to be changed to accommodate less 

demanding design conditions. Usually this can be 
performed only step by step via an iterative procedure. The 
steps can contain model simulation and also experiments 
on the real process. Therefore the optimization is often 
termed the redesign of the reference model. In the case of 
low-order reference models the time constant of the model 
(i.e., the bandwidth) can be determined by explicit design 
expression if the process model and the amplitude 
constraint U are known. 

max 

Let the process be given by the pulse transfer function 
0.125z‘‘(l + 1.6z -1 ) 

G = 1 - fV? -IT (28) 

l-0.5z -1 j^l-0.8z -1 j 



Apply a combined iterative identification and control 
test [12]. The following reference models of unity gain are 
used for the design 


K pt = arg 



= arg 


\ min 

i *» 



0.9z~ 
1-0. lz _1 


and R v 


0.9z -1 
1 - 0. lz _1 


(29) 


where the chosen criteria 7Ls = |P _ T-| an d 

Ties = ||l - T, | state that each reference model should 

approach the ideal unity. This task must be solved under 
the constraints for the regulator output u dJ . Here U 

means the allowable region for u , i.e., 1/: |w| < 1 (see 
(24)). 

Redesign of the Reference Model 

The optimization task (27) is very difficult because the 
solution is always on the border of the limited region. 
There is no analytical solution except for some low-order 


At the start of the iteration the model 
0.1z _1 (i + 4.0z _1 ) 

G ° = 7-W- i\ ( 3 °) 

(l-0.2z _1 |(l-0.9z _1 ) 

is assumed. A square signal with periodic time of 40 
samples is applied as reference signal. In the simulation it 
is assumed that the additive noise y n is white noise, 

whose variance is a = 0.01 . The number of the 

U 

processed samples is A = 100 . Because of the small 
output noise the identification is performed by a simple 
off-line LS method. The regulator is designed by the YP 
method, assuming an IU process. 
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actuation by the end of the iteration. The obtained 
redesigned reference model is 


- 0.5022z~ 1 

r 1 - 0.4978z _1 


(31) 


Figure 13 shows the time function of the output of the 
reference model (continuous line) and the closed system 
(dashed line). Fig. 14 presents the BODE diagram of the 
original (continuous line) and the redesigned (dashed line) 
reference model. The prescribed over-excitation could be 
reached only by substantially slowing down the closed 
system. 



Figure 14. Bode diagram of the original and the new reference model 


Figure 12. Response of the YP regulator after the iteration 

The output of the regulator is presented in Fig. 11, where 
it is seen that the over-excitation is very high at 900 %, i.e., 
u t = 9 . Assume that the actuator can realize only u t = 5 . 

This requires the redesign of the reference model R r . 



Figure 13. Output of the reference model and the process 

The condition u t < u t = 5 needs to be prescribed for the 

reference model redesign iteration. It can be seen in Fig. 12 
that the control output is according to the prescribed over 


The final conclusion is that the control quality can never 
be prejudice free, more exactly never process independent 
and never the installed actuator independent property of the 
final control system. 



Figure 15. Illustration of the inequality of (9.5.1) 

4. The Heisenberg-T y pe Uncertainty 
of Control 

The condition of robust stability (18) already contains a 
product inequality. Here |7 T (y<i>)| (although it is usually 
called a design factor) can be considered as the quality 
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factor of the control. The other factor, however, can be 
considered as the relative correctness of the applied model. 
In the light of practical experience control engineers favor 
applying a mostly heuristic expression 

(quality of the control) x (robustness of the control) < limit 

This product inequality can be simply demonstrated by 
the integral criteria of classical control engineering. Let I 2 

be a square integral criterion (Integral Square of Error: ISE ) 

* 

whose optimum is I 2 when the regulator is properly set, 
and the Nyquist stability limit (i.e., robustness measure) is 
p m . The well-known empirical, heuristics formula is 

I* 

j- P m ^ limit (32) 

A 

The inequality is illustrated in Fig. 15. The fact that the 
quality of the identification (which is the inverse of the 
model correctness) can have a certain relationship with the 
robustness of the control, is not very trivial. This can be 
observed only in a special case, namely in the identification 
technique based on KB parameterization [1] [2] when 
l&jp = -#< . Introduce a new relationship for the 
characterization of the quality of the control 

' ,33) 


)-\ 

-e(yco) 


1 


1 

J" 

>’n (i®) 


1 + CP 


i+Z 


Notice that 5 is the absolute value of the sensitivity 
function. Obviously, 8 p = 1 for all frequencies (here 

p = |l + Z| ). Of course, the same equalities are valid for the 
minimum and maximum values, i.e.. 


Pm(c') = min p(co,c) =min(p) 

' ' CD L ' > A 00 


5 M (c) = m ax 5 

\ ’ rn 


(»4 ! 


; i.e., their 


product is unity: S M P m = 1 (34) 

Denote the worst value of these measures by 


Pm = m 4X 

min 

p(co,c) 

> = max min(p) 

c 

, CD 


C L CD 

- m i n I 

max 

s(co,c) 

| = min max (8) 

c 1 

k CD 

V /_ 

j C L to 


Pm 

The above three basic relationships can be summarized 
in the inequalities below 

8p = 1;8 M p m = 1; 8 M p m =1 (36) 

where the following simple calculations prove the 
existence of (34) and (35) 


p m (C) = min 1 + CP\ - - 


1 


1 


1 + CP 


1 + CP 

fT 

m 

8 


u(c ) 


5 m = min < max 


c co 1 + CP 


max < max-— 

c ® 1 + CP 


-\j max |min|l + CP|j = l/p m (38) 

Given (34), (35) and (36) further basic, almost trivial, 
inequalities can also be simply formulated 


> Pm(^)-Pn 

^ = 8 m ^8 m (c) = -2_ 


5 M (c)“ FmV 

The above results are not surprising. The fact, that they 
are valid even for the modeling error in the case of 
^-parameterized identification methods makes them 
special. So it can be clearly seen that when the modeling 
error decreases, the robustness of the control increases. 

Namely, if the minimum of the modeling error 8 M is 
decreased, then the maximum of the minimum robustness 
measure p m is increased, since 8 M p m = 1. 

Similar relationships can be obtained if the J~C 2 norm of 
the “joint” modeling and control error is used instead of the 
absolute values. Introduce the following relative fidelity 
measure 

||£tta ||_ Ill'll/-* II II 

° = n = n ; W (40) 

IWI 2 IWI 2 

The upper limit for this measure can be formulated as 

CT = Tnr = FT _ “ = 5 m(^) (4i) 

Ik n II 2 Ik n II 2 ^ + CP 

so it is very easy to find similar equations for a . Let 


a M = min 

c 


= max a(/,c) and 
linjmax j . 


Using these definitions and the former equations we 
obtain the following interesting relationship 


a M - a M 


(C)S8 M (C) = 


for the relative quadratic identification error. 
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Use the first-order reference model 


R 


n 


b nl Z 


-1 


l + a nl z 


-1 


(43) 


for the design of the noise rejection in the IS process. Here 
the maximum of the robustness measure is 


Pm=Pm,is = 0 - 9 according to 
Pm, IS = Pm, IS p = 0) = |T - 


= min 

CO 


1 


1 -IL 


(44) 


and 


Pm, IS _ 



(45) 


This gives a robust measure of p^ IS =0.9 under the 
parameters b nl = 0.2 and a x = -0.8 chosen for 
reference models with unity gain (same as in (29)). 

The values of the typical variables (see above) are 

^'iF'irhr 1 - 1111 * 8 «P" =1 (46) 

Pm rm l \ 


a 


o _ 
M ~ 



/ 2 =1.054 => 

Vki-i| 


®° M Pm= 0.9486 <l° M p° m =l (47) 



Figure 16. Illustration of uncertainty relationships (41) 

Considering the data of (46) and applying again the 
relative sampling time x = T s /T n , the different measures 

in (42) are illustrated in Fig. 16. Here T n is the time 

constant of the continuous time (CT) first-order reference 
model. 


Introduce the following coefficient for the excitation 
caused by the reference signal 


5= 


hi . 

• M 

y„ 2 



*0 


(48) 


which represents a signal/noise ratio. Investigate the 
product g p (which is called the uncertainty product) in 

an iterative procedure where the relative error [E of the 
model is improved gradually. For simplicity, let us assume 
an IS process. It can be simply derived that 


Pm ^ Pm ^ 


1+i WI 


:a o Pn 


minll+7^1 l-^||^ n || 
= 1 


^->o 


(49) 


i.e., 


ap m < ap° < a 0 Pm - 1 or 
^-»o 


a ^ = 7^- = 8 m ^ 8 M = - 

Pm p 


(50) 


g 0 = <j(£ = 0) . Similarly to the notations 
and g m applied above, the notations 


where 

ct m(c) 

a m (f) = min <j^£,C^j and g^ = o m (£ = 0) can also be 
introduced. It is not an easy task, however, to derive the 
relationship between g^ and g q or g m and G M ^cj. 

The simplest case to investigate (50) is when £ = 0, since 
then 


< (c)<5^(c)- (51) 

This equation gives a new uncertainty relationship, 
according to which 

!f=-min |l + Cp\ < 1 (52) 

c I U ->0 

\\2 

The product of the modeling accuracy and the robustness 
measure of the control must not be greater than one, when 
the optimality condition £ - 0 is reached. The obtained 
uncertainty relation can be written in another form, since 

min|l + CpI =1 (53) 

c I \t=o 



The earlier results of control engineering referred only 
for the statement that the quality of the control cannot be 
improved, only at the expense of the robustness, so this 
result, which connects the quality of the identification and 
the robustness of the control, can be considered novel. 

This phenomenon can arguably be considered as the 
HElSENBERG-type uncertainty relationship of control 
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engineering, according to which 


Here Az and Ap are the alterations of the canonical 
coordinate and the impulse variables, respectively, and 
thus their inverse corresponds to the generalized accuracy 
and “rigidity” which are known as performance and 
robustness in control engineering. 

The consequence of the new uncertainty relation is very 
simple: ^-parameterized identification is the only method 
where the improvement of the modeling error also 
increases the robustness of the control. With other methods, 
and other identification topology, modeling and control 
errors are interrelated in a very complex way, and in many 
cases this relation cannot be given in an explicit form. This 
is the main reason why it is difficult to elaborate a method 
which guarantees, or at least forces, similar behavior by the 
two errors, though some results can be found in [3]. 



Figure 17. Relationship between the control and identification error in 
the general case 



Figure 18. Relationship between the control and identification error in 
the case of the ^-parameterized identification method 

There is a myth in the literature concerning the 


antagonistic conflict between control and identification. A 
“good” regulator minimizes the internal signal changes in 
the closed loop and therefore most of the identification 
methods, which use these inner signals provide worse 
modeling error, if the regulator is better. The exciting 
signal of ^-parameterized identification is an outer signal 
and therefore the phenomenon does not exist. The relevant 
feature of this relationship is shown in Figs. 17 and 18 for a 
general identification method and a ^-parameterized 
technique. 

In Fig. 17, there is no clear relation between 5 ID and 8 , 
or g id and a, and therefore there it is not guaranteed 
that minimizing 8 M increases p m . In Fig. 18 5 ID = 5 
and a ID = a , and thus the minimization of 8 M directly 
maximizes p m . Thus if during the iterative identification 

the condition \\i J| = 0 is guaranteed then, at the same 

11 * 1,00 *->00 

time, the convergences 5^ = 5^ and p^ = p^ are 

»oo »co 

ensured. 


5. Irregularities in Classical Control 
Methods 

A further myth in the control literature is that everything 
is right and errorless in the classical works of theory. This 
is unfortunately (or fortunately ?) is not right. T. Keviczky 
recognized that the solution of the classical LQR paradigm 
does not provide full pole placement and some areas can be 
considered as irregularities of this classical theory. 
Keviczky and Banyasz [9], [10] gave a detailed analysis 
proving that there are some poles, it is interesting that the 
slower poles, which can not be allocated by the classical 
methodology. In some earlier references it was also studied 
that some poles are unreachable in the LQR theory by 
Johnson [11]. Finally Bokor and Keviczky [12], [13] 
presented a possible method, which solves the irregularity 
of this paradigm. 

The LQR (Linear system - Quadratic criterion - 
Regulator) problem 

This optimization paradigm was formulated by the 
general, quadratic criterion [7], [8] 

I = T\W{t)W x x{t) + w n u 2 {t)Ut (55) 
Z 0 

where is the state vector, u(t^ is the input of the 

process, respectively. The positive definite W x stands for 

penalizing the variations in the state space, w u is for 

penalizing the energy of the control action. The classical 
solution, minimizing (55) is a negative state feedback (, SF ) 
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'(0=-*LQ*(0 


where k is given by 


kL - -t T p 

w 


<56) 

the equation (63) can be rearranged into a new form 

(57) M(^)]R(-5) = [a(5) + A: t 'P(5)6][a(-5) + A: t ^(-5)6] = 


The symmetric positive semi definite matrix P can be 
obtained from the solution of the algebraic RlCCATl 
equation [4] 


PA + A T P-—Pbb T P = -W (58) 

W xv/ 

U 

Analytic solution is not possible, because this equation is 
nonlinear in P , therefore only numeric solution can be 
obtained by MATLAB and other CACSD programs. 
Introducing the orthogonal factorization 

W K = G t G (59) 


= AES’ 


-s)+[g T 'F(s)z>][gN'(-s)z>] 


( 66 ) 


which provides the quadratic polynomial solution of the 
KALMAN equation. Thus the final quadratic equation, 
ensuring relationship between the process A ( 5 ), design 

M(vS') and weighting G(s) polynomials, is 

= A ( 5 ) A (- 5 ) + G(s)G(-s) 
or |m(s )| 2 =|a(s)| 2 + |g(s)| 2 


the closed-loop system is stable if the auxiliary process 

v = Gx (60) 

is observable. 

The Frequency Domain Solution of the LQR Problem 
The LQR approach is widely used for control problems 
in all over the world, however, in a practical problem it is 
not an easy task to find the best W x and w u weights, 

which are usually obtained by trial and error iterative 
methods. The LQR problem has an almost forgotten 
frequency domain solution, too, which will give us a 
deterministic design process to find useful relationships 
between the classical pole placement SF solution and the 
LQR paradigm. It can be shown that the simpler dyadic 
factorization [7] 


or in the general form 

w u |m(s)| 2 = % |A(s)| 2 +|G(s)| 2 (68) 

Observe that the solution tends to M(s) = A(s) if 

w u —» 00 and g T x = 0 if w u — » 0 . Do not forget that 

K(s) and G(s) are (n - l)th order [8]. 

Some Anomalies in the LQR Problem 
The solution of the polynomial equation can be a direct 
coefficient comparison or a spectral factorization approach 
[8]. Consider some examples in the sequel. 

Example 1 

Consider a first order example with 


W x =gg j \ g=[_g v ...,g n _ vgn J (61) 


can also be used. The frequency domain condition of the 
minimum of (55) is called the KALMAN equation [7], [8] or 
sometimes it is named frequency domain identity 


i+& T <&(s)& 


= w + 


g 


®(s)L 


(62) 


Assuming unity weight w u = 1 the equation becomes 
even simpler 

|i+* t ®(s)6| 2 =i+||^ t ®(5)6|| 2 =i+|s t ®(s)6| 2 (63) 

Using the well known relationship of complex functions 

{Z(s)} 2 = \z(sf =Z(s)Z(-s) (64) 

and introducing the (n - l)th order polynomial G(s) as 
the numerator of 


A(s) = s + a l ; M(^) = ^ + r 1 thus G (s) = g x (69) 

The two sides of (67) are 

-s 2 - r x s + r x s + r 2 = - s 2 - a x s + a x s + a 2 + g 2 (70) 
and the solution is 

if = a; +g 1 ;>o (7i) 

and 

k x =r x -a x = + gf - a x > 0 (72) 

If we want to ensure (place) a required pole then the 
necessary weight in the LQR problem is 

g, = Jrf-af (73) 

It is easy to see that only such q can be placed, which 
fulfills the condition 
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2 2 


K > \(L 


> r x > \aA 


(74) 


for stable design polynomial M(v). 

So this example shows that only a faster pole can be 
placed by the LQR optimization comparing to the original 
process pole. 

Example 2 

Consider a second order example with 

A(v) = s 2 + a x s + a 2 ; M(v) = s 2 + r x s + r 2 
thus — g x s + g 2 

The two sides of (67) are now 


(75) 


( 5 


+ r x s + r 2 j - r x s + r 2 j = + a x s + a 7 

p-a 1 ^ + < 3 2 ) + (g 1 ^ + g 2 )(-g 1 ^ + g 2 ) 
and the solutions are 

r 2 — ^ a 2 <§2 ^ a 2 
r \ = 2 ) + ( a l 2 +^ 1 2 ) = 

= ^2(r 2 -a 2 ) + (a 1 2 +gi 2 ) > flj > 0 

The £F to be applied is given by 

^ > 0 ; k 2 = r 2 - a 2 > 0 

For pole placement the necessary LQR weights are 


» 


§2 


/r 2 2 - a\ 


(76) 

(77) 

(78) 

(79) 

(80) 


and 


gi Jr 2 2^r 2 J2 (p r (81) 


where 


2 r l ~ 2r 2 
H'r ~ 




(82) 


and 


2 a i 2a 7 
Ma =- 


(<) 2+ U 


(83) 

2 2 
It is easy to see that there are such |r p r 2 | domains, 

which can not be reached by any -jgj, g 2 | selection !!! 
These conditions are 


2 2 
r 2 > a 2 


\r n > \a~ 


• > al 


(84) 



(85) 



Figure 19. Unreachable design parameter domains 

These conditions are graphically demonstrated on Fig. 19, 
where the shaded area shows the unreachable design 
parameters for the case of open-loop process parameters 

|^ 2 1 = 0.8 and 2p^=0.5. 

One can check these results either via the solution of the 
RlCCATl equation (very time consuming method) or by the 
spectral factorization approach 

M(s)M(-s) = [a(s)A(-.s) + G(-s)G(.s)] + 

r ( 86 ) 

[a (s) A (- 5 )+ G(—s)G(.S')] 

as the solution of (67), i.e., by 

R(v) = [A (s) A(-v) + G(v)G(-s)] (87) 

Solutions for LQR-pole placement 
A new criterion 

7= 2 J u T (*K* ( f T* T ( ? Ku m 0+ w y (Or < 88 ) 

was proposed by Bokor and Keviczky, which solves the 
problem and offers a possible way to overcome this 
anomaly of the standard method. 

The full pole-placement solution can be obtained only if 
the criterion (86) is used, which penalizes interaction of the 

state vector and the control action by the x T (H4) 
cross-term. The optimal state feedback can be provided by 

k = -w^ 1 (w m + P b) or k T = -wT 1 (w ux + * T i>) (89) 


and 


introduced in [11] and [12]. 
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6. Conclusions 

The purpose of this paper is to highlight some interesting, 
may be philosophical, paradigm of modeling and control as 
the extended abstract indicated. 

Such special and not widely known problems are 
discussed here, which are worth further study and 
investigation. 

“I believe that the progress of science should be rather 
measured by the raised and not by the solved problems !” 
as Eddington stated !!! 
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